Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 10 April 2013 (MN KTeX style file v2.2) 

Stokes imaging of AM Her systems using 3D inhomogeneous 
models - II. Modelling X-ray and optical data of CP Tucanae.^ 



10 April 2013 



^ ■ K. M. G. SilvaH, C. V. Rodrigues\ J. E. R. Costa\ C. A. de Souza^J, D. Cieslinski^ 
O & G. R. Hickel2§ 

^ ^ ' ^ Instituto Nacional de Pesquisas Espaciais/MCT — Av. dos Astronautas, 1758 — 12227-010 - Sao Jose dos Campos - SP - Brazil 
5-H ^ Universidade do Vale do Paraiba-Av. Shishima Hijumi, 2911 - Urbanova - Sao Jose dos Campos - SP — Brazil- 

6 



> 

in 

o 
\o 

(N 

^' 

O 
en 



X 



ABSTRACT 

The viewing geometry of the polar CP Tuc that better explains its optical and X-ray 
light curves is controversial. Previous modelling of white-light polarimetric data considered 
the partial self-eclipse of an extended inhomogeneous emitting region. Alternatively, phase- 
dependent absorption has been used to reproduce the X-ray data. This paper presents new 
optical polarimetric data of CP Tuc and a model that consistently explains its optical and 
X-ray data. The model was based on an extension of the CYCLOPS code that added X-ray 
bremsstrahlung emission and pre-shock region absorption to the original version, which only 
accounted for cyclotron emission. The new code creates the possibility of simultaneous optical 
and X-ray fitting. We show that self-eclipse and absorption data have distinct signatures on 
the X-ray spectra. Although we were able to reasonably fit the CP Tuc optical data to cases of 
absorption and self-eclipse, we were only able to reproduce the X-ray orbital modulation after 
considering the absorption in the pre-shock region. Specifically, we were unable to reproduce 
the X-ray observations in the self-eclipse case. We found that the primary emitting region in 
CP Tuc is located near the rotation pole that approximately points to the observer 
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1 INTRODUCTION 

Cataclysmic variables (CVs) are compact binaries in which a low- 
mass late type main-sequence star transfers matter to a white dwarf 
(WD) due to Roche lobe overflow. Magnetic CVs (mCVs) are sys- 
tems in which the WD possesses a strong magnetic field that con- 
trols the dynamics of the accretion flow near the WD. mC Vs are di- 
vided into two subclasses: AM Her systems (i.e., polars) ( Croppei 
19901) and DQ Her systems (i.e., intermediate polars) JPatterson 
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1994) . The typical magnetic field on the WD surface of polars and 



intermediate polars ranges from 5-100 MG andO. 1-lOMG, respec- 
tively. 

The WD spin of polars is synchronised with the orbital period, 
and the magnetic field prevents the formation of an accretion disc 
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( ICroDpeiJll990l) . The magnetic pressure overcomes the accretion 
flow pressure at the Alfven radius. In this so-called coupling region, 
the flow from the secondary star leaves the ballistic trajectory and 
turns to a magnetically driven stream to the WD surface, thereby 
forming an accretion column. Near the WD, the material reaches 
a supersonic velocity, and a shock front is formed. The region be- 
tween t he shock front and the W D surface is called the post-shock 
region (Lamb & Masters 1979). As a result of the shock, this re- 
gion is hot, with temperatures greater than 10 keV. The region cools 
down via cyclotron and bremsstrahlung emissions, which dominate 
the optical and X-ray ranges, respectively. 

The cyclotron emission from the post-shock region accounts 
for the orbital variation of the optical flux and the polarisation 
of polars. Recently, we develope d a code to re produce the opti- 
cal emission of mCVs: CYCLOPS JCosta & Rodr igues 2009). This 
code adopts a 3D-representation of the accretion column and WD, 
which allows us to study the complex physics and geometry of 
mCVs. This code calculates the cyclotron emission from inhomo- 
geneous post-shock regions and the attenuation from the pre-shock 
region by considering the Thomson scattering. .Costa & RodriguesI 
( 120090 briefly review previous models. 

The X-ray emission of polars is consistent with hot plasma 
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emission i Lamb & Mastersll 19790 . The Raymond-Smitii and Meka 
models JRavmona 1 20091) can reproduce this emission. Initially, 
the X-ray data of mCVs were modelled using a single tempera- 
ture hot plasma continuum and Gaussian profil es to represent the 
few iron emission line s observed in the spe ctra (llshida et alj|l99ll : 
lEzuka & Ishidi 1 19991 : iTerada et al.l 1200 ih . However, the X-ray 
spectra of certain mC Vs cannot be fit with a single bremsstrahlung 
component, which indicates the presence of mu lti-temperature 
plasma in the post-shock region dEzuka & Ishidalll999) . The ra- 
diative co oling in the post-shock region gene r ates the temperature 
struc ture CWu, Chanmug am & Shavivl Il995t iDone & Magdziara 
[T998 : Cropper et al. 1999). 

Absorption is necess ary to reprodu ce the spectral energy de- 
pendence of mCVs (e.g.. lMukaill2011r) . On one hand, interstel- 
lar absorption in the 0.1-2 keV region with typical column den- 
sities of_J_xJ.O^i£Sl~f_is observed in the X-ray mCV spectra 
(e.g. lEzuka & Ishidal 1 19991) . Moreover, an additional absorption 
varia ble in the orbital cycle has been observed in certain s ys- 
tems JDone & Magdziarzll 19981 : ICropper. Wu & Ramsavll2000l) . A 
partial occultation of the emitting region by the neutral material 
in the pre-shock region causes this absorption. Absorption mod- 
els adopt an ad-hoc representation of these geometric al effects 
JDone & Magdziarzlll998l : ICroDper. Wu & Ramsavll2000l) . 

In addition to absorption, the occultation of the emitting 
region behind the WD li mb can produce orbi tal variations in 
the emission of m CVs (King & Shavivj 1 1984) . For instance, 
lAUan. Hellier & Beardmore (1998) discussed these two processes 
in the context of the phase-dependent effects observed in the inter- 
medi ate polar EX Hya. B oth effects are possible in AM Her sys- 
tems dWatson et alj 19891) . The absorption by the accretion column 
can produce narrow dips due to the material located near the cou- 
pling region, broad dips due to the material near the post-shock re- 
gion, or both dWarren, Sirk & Vallergalll995USirk & Howelllll998l : 
iTovmassian et al.ll2000l) . 

The X-ray fits of polars are usually conducted without consid- 
ering the optical data, partially because distinct codes are used to 
model each s pectral range. One e xception is the polar RX J2115- 
5840 model dRamsav et al.ll2000l) . In this case, the cyclotron 2D- 
emitting region, obtained from the optical data fit using the method- 
ology of Potter, Hakala & Cropper ( 1998), is used to define the lo- 
cation of the X-ray emission. However, the density and temperature 
profiles adopted when fitting each energy band are not identical. 

CP Tuc is a polar that presents energy-d ependent orbital X- 
ray m odulations in the high state of brightness dMisaki et al.ll 19951 
1 1996). The ASCA satellite discovered CP Tuc in February 1995 
dMisaki et al.lll995l) . Two suggestions have b een made regardin 
the nature of this modulation (as for EX Hya). iMisaki et al.l dl99i 
fit the phase-resolved s pectra using a part ial-covering absorption 
model. Alternatively, iRamsav et al.l d 19990 argued that accretion 
stream absorption was an unusual mechanism to modulate the X- 
ray light curves of polars and considered a self-eclipse scenario to 
be more likely. These authors fit the optical polarimetry data using 
one self-eclipsed region and suggested that a gradually eclipsed, ex- 
tended inhomogeneous region can explain qualitatively X-ray mod- 
ulation. 

This paper presents new optical polarimetric data from CP Tuc 
and extends CYCLOPS to model optical and X-ray data together. 
This new version of CYCLOPS is then used to study the spectral 
signature of the X-ray modulation via absorption or self-eclipse 
as well as how these situations can be distinguished using high- 
energy observations of polars. Subsequently, CYCLOPS was specif- 
ically applied to CP Tuc. This work is organised as follows: Sec- 



tion |2] details the modifications conducted on CYCLOPS. Section|3] 
presents a study on the phase-dependent effects in the X-ray spec- 
tra of AM Her systems caused by self-eclipse or partial covering by 
the pre-shock region. Section |4] presents the model of the CP Tuc 
optical and X-ray observations, which includes new optical polari- 
metric data. The results are discussed in light of the previous litera- 
ture. Section|5]presents the conclusions. The preliminary resul ts of 
this study were present ed in lSilva. Rodrigues & Costal d20Ilah and 
lRodriguesetal.ld20Ilh . 



2 EXTENDING CYCLOPS TO X-RAYS 

CYCLOPS is a code designed to reproduce the cyclotron emission 
that originates in the post-shock regio n of mCVs. A deta i led de - 
scription of CYCLOPS is presented in ICosta & RodriguesI d2009h . 
This section describes the implementation of the X-ray emission 
in CYCLOPS. Initially, we present an overview of the code (Sec- 
tion 12. It . Subsequently, we outline the implementation of the 
bremsstrahlung emission (Section [2. 2) . the absorption that occurs 
in the accretion column (Section|23}, and the procedure for simul- 
taneously fitting optical and X-ray data (Section|2!4j- 



2.1 CYCLOPS 

CYCLOPS builds a 3D grid to represent the entire accretion column 
from the threading region to the white dwarf surface. The entire 
column is divided into several volume elements (voxels) accord- 
ing to the chosen spatial resolution. The lines of a centred dipolar 
magnetic field define the geometry of the accretion column. The 
magnetic axis has an arbitrary direction relative to the white dwarf 
rotation axis. Two footprints (northern and southern) are on the WD 
surface. This paper considers only one of these emitting regions. 
We defined the post-shock region from the WD surface to the height 
of the shock front. The emission comes from this region. The pre- 
shock region extends from the shock front to the threading region. 
The phase-dependent absorption is produced in the pre-shock re- 
gion (see Section [231 l. The functions described below define the 
electron density and temperature of each voxel of the post-shock 
region. The magnetic field strength in each voxel was calculated 
using the magnetic dipole. 

A Cartesian frame was defined for each orbital phase in which 
the emission was calculated. The z-direction of this frame was the 
observers line of sight. Therefore, the observer views the accretion 
column as a 2D set of sight lines. The radiative transfer for each 
line of sight was solved by considering the emission and absorption 
processes that were implemented in the code. 

The CYCLOPS input parameters are as follows: 

• j is the inclination of the orbital plane of the system relative to 
the observer. 

• /3 is the angle between the rotation axis and the centre of the 
northern region. 

• /; defines the tangential location of the maxima of the density 
and temperature. It ranges from to 1. /; = 0.5 stands for a central 
location. 

• A_B and Along define the size of the coupling region and, 
therefore, the footprint of the emitting region on the WD surface. 
The coupling region is defined as a 2D region that contains a ref- 
erence point located at a radius Rth from the WD surface and at 
longitude longth- Rth and longth are defined by /? and the mag- 
netic field geometry, see description of the geometry of the emitting 
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region in ICosta & RodriguesI 1 120091) . The coupling region radial 
length is 1 ± Ah Rth, therefore A/j varies from to 1. Its lon- 
gitudinal extension is longth ± 2A!on9 (1-/0. ^long varies from 
to 180 degrees. 

• /i is the height of emitting region in WD radius units. 

• Bpoie is the intensity of the magnetic field on the magnetic 
pole. 

• Biat and Biong define the direction of the dipole axis. 

• Tmax and Nmax ate the maxima of the electron temperature 
and density, respectively. 

We chose arbitrary functions to describe the variations in tem- 
perature and density with regard to height of a specific voxel from 
the WD surface, h^ox, which are 

l^vox 

T 



J- yil'vox ) — J- max 6Xp 



(2.5 



\ h J 
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N(hvox 



N„ 



:exp 



-2.5 



f^vox 
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These functions resemble the density and temperature struc- 
tures in the h ot and magnetised post-shock regions obtained by 
ICropper et all ([7999) (see their Figure 1) and Saxton et al. (2007) 
(see their Figure 2). However, Equations [T] and |2] do not represent 
shock solutions and are not necessarily consistent with other model 
parameters. 

With regard to the density and temperature profiles perpendic- 
ular to the radial direction, the code allows us to consider a con- 
stant profile or modify T[hvox) and N{hvox) using the following 
expressions: 



T{d) = T{h^ox) e 
and 

N{d) = N{Kox) e 



-Vd 



-Vd 



(3) 



(4) 



where d is the distance from the voxel to the reference point in 
which the temperature and density values are greatest. These ex- 
pressions are plausible representations of the tangential decay of 
temperature and density in the accretion column. They are useful 
for testing whether the tangential variation of density and tempera- 
ture affects the mCV emissions. 

The first version of CYCLOPS only considered the cyclotron 
emission and bremsstrahlung absorption processes; moreover, it 
was used to reproduce the optical and infrared polarisation of 



was usea to reproauce tne optical ai 
V834 Cen dCosta & Rodriguesi2009l) . 



2.2 Bremsstrahlung emission inCYCLOPS 

The bremsstrahlung emissivity of a fully ionised magnetised hy- 
drogen plas ma (A^e = Nio„) was calcul ated for each voxel (k) 
according to lGronenschild & Mewd ( 1197 81) : 



JE,k = 1.032 X 10"" g nI E~^ T" 



(5) 



where A'^e {cm~ ) is the electron number density, E is energy, 
T(K) is temperature, and g is the n on-relativistic Gaunt factor 
(IMewe. L emen & v an den Ooro 1986). 

The emission of a given orbital phase is the sum of the fluxes 
from all optical paths (see Section l2Tt . The radiative transport is 
calculated for each optical path from the farthest to nearest vox- 
els by considering the incident radiation from the previous vox- 
els. Physical quantities vary along the optical path following the 



equations for density and temperature presented above. This vari- 
ation naturally incorporates the shock structure and optical depth 
effects, which are im portant to correctly calculate the emerging 
flux (see lCroDper et al. 1999). These opacity effects in X-ray wave- 
lengths are important when the column density is greater than 
1.3 lO^'^ cm~^. The X-ray models use at least 64 voxels in the ra- 
dial direction to avoid under-sampling the temperature and density 
functions. 

Considering the thermal case, which is appropriated to mCVs, 
the specific intensity for a given energy E from each line of sight / 
at orbital phase ph is 



lE,i,ph = 2_^(^E,Lph,(k~i) - bbE,ph,k) e ^'^ + bbE,ph,k , (6) 
fc=i 

where bb is the Plank function and nk is the number of voxels k at 
line of sight I. Tk is the optical depth of a specific voxel, which is 
calculated as 

JE,k ,„, 

Tfc = 7T S, (7) 

bbE.k 

where s is the optical path equal to the voxel length. lE.i,ph can 
be used to construct the image of the post-shock region in a given 
orbital phase using the frequency selected. In the absence of ab- 
sorption, the total observed flux at the orbital phase ph is 

2_^^E.l,ph, (8) 



FE,ph = C * y^ lE,l.ph X 

I I 

where (7 is a constant, A^ is the emitting area, and D is the distance 
to the system. 

2.3 Absorption 

CYCLOPS can account for two sources of absorption: from the in- 
terstellar medium (ISM) and from the pre-shock region. 

To calculate the photo-absorption cross se ction for both situ - 
ations, we used the IDL routine bamabs (Kashv ap & Drakdl 20001 ). 
This absorption cross section is a function of temperature and abun- 
dances. We considered a homogeneous material with solar abun- 
dance and temperature equal to 10000 K. 

The interstellar absorption is constant along the orbital cycle 
and depends on the quantity of the interstellar material in the direc- 
tion of the source. The interstellar column density, Nl^i''' (cm~^), 
was estimated using FTOOL Nh Q NI^/'' was used as the upper 
limit of the interstellar column density of hydrogen in the fit. 

The pre-shock region followed the magnetic field lines from 
the post-shock to the threading regions. All the voxels in the pre- 
shock region were assumed to have the same density and an upper 
limit set by the minimum density of the post-shock region (see be- 
low). The temperature was also constant. 

The pre-shock material absorption is phase dependent. This 
absorption was considered for the phases in which the pre-shock 
region fell between the observer and the post-shock region. CY- 
CLOPS included the absorption for each line of sight I, which were 
summed to produce the total flux coming from the object in a given 
orbital phase i. In particular, each line of sight has an absorption 
rate that is proportional to the column density of the pre-shock re- 
gion. 

Considering the absorption from the ISM and the upper por- 
tion of the column, the flux observed in each orbital phase i is 

^ |http: //heasarc .gsfc . nasa.gov/cgi-bin/Tools/w3nh/w3nh .pi) 
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E.l,ph e 
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(9) 



where Tpr^ = h s cjph N^^" , tmi = cTph iVf,,^", (Jph is the 
photo-absorption cross section, N^^" is the minimum number den- 
sity in the post-shock region, and ki is the number of voxels in the 
pre-shock region for each I. The constants A and B are values be- 
tween and 1 that defines the fractions of N^*" and N'Ji'^' , which 
are used to fit the data. 



2.4 Fitting data 

Our aim is to simultaneously fit the X-ray and optical data of 
mCVs. To determine the model parameters that best fit the obser- 
vational dataset, we used the least value of the figure of merit, y^, 
which was defined as 



X = xlpt + addxlx, 



(10) 



where Xopt measures the agreement between a given parameter set 
and the optical data, and Xrx plays the same role for high-energy 
data. We used light and polarisation curves in the optical regime, 
and we used spectra for the X-ray energies. The factor add allowed 
us to choose the relative contributions of the optical or X-ray data 
in the fitting procedure. This flexibility is important for achieving a 
satisfactory fit in both energy ranges. 

The code automatically defined the wavelengths and the ph 
orbital phases for which the model should be calculated using the 
input datasets. Usually, the orbital resolution of optical data is supe- 
rior than that of the X-ray spectra. Consequently, a X-ray spectrum 
was produced for all ph orbital phases of the optical dataset. These 
spectra were summed over specific orbital phase intervals to obtain 
the spectra in the j interval phases that correspond to the observed 
X-ray spectra. 

The figures of merit slightly differ in the two spectral domains. 
Specifically, we calculated the functions xl;„ j and Xb^x using: 



2 
XEopt 



and 



2 

XEr 



Eopt ph 



EEi' 



,ph 
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(11) 



(12) 



where Eopt and Erx are the energy in optical and X-rays wave- 
lengths, respectively, Fe^^,j and FE„pt,ph are the fluxes calcu- 
lated using the core of the code, and dE^^j and dE„pt,ph are the 
observed (total or polarised) fluxes. The multiplicative value fcyc 
was used to normalise the model to compare with the data. This 
procedure was conducted because the distance and the size of the 
emitting region were not used to correctly scale the model. We esti- 
mated fcyc using the optical data. The additive constant F"*" repre- 
sents the non-polarised emission in optical wavelengths. This emis- 
sion originates from the stellar components of the system. Each fre- 
quency has a different F"^^ value because this component is energy 
dependent. The superscript conv in the X-ray flux is explained be- 
low. 

The X-ray flux, as calculated by CYCLOPS, Fe^^^ph, cannot 
be directly compared with the data; we must consider the response 
of the detection system to fit X-ray data. During the data reduc- 
tion, we obtained two files: the detector redistribution matrix file 
(RMF), which accounts for detector gain and energy resolution. 



and the ancillary response file (ARF), which accounts for the ef- 
fective area of the telescope/coUimator (including vignetting), fil- 
ter transmission, detector window transmission, detector efficiency 
and any additional energy dependent effects (see OGIP Calibration 
Memo CAL/GEN/92-OO Jfl). We used the IDL routine cojivrmf 
JKashvap & Drakdl200Ci) to convolve the original flux from CY- 
CLOPS, FEr^,ph, with the response matrices to obtain the flux for 
observation comparisons, Fe°^ 



,ph • 



F'Conv 
E. 



ph 



{ARF FE,,,ph) X RMF. 



(13) 



In addition, three other internal parameters were automatically 
optimised during the fitting procedure: 

• an offset in the model orbital phase, Sphase', 

• the fraction of the upper limit of the electron density in the 
attenuation region, A (see Section lZJl l: 

• the fraction of the upper limit of the extinction due to ISM, B 
(see Section [23] l. 



As in the original code JCosta & Rodriguesll2009l) . the func- 
tion x^ was mini mised using two alg orithms. First, we used the 



pikaia algorithm dCharbonneaul 1 1 995h with a broad range of pa- 
rameters. Then, the best solution s were refined using the amoeba 
algorithm (e.g.. |Press et al.lll992h . 



3 ORBITAL VARIATION OF THE X-RAY SPECTRA 

The orbital variation of high-energy spectra is a common feature in 
polars (see the Introduction). CP Tuc presents this type of variation 
and our primary aim was to reproduce its X-ray spectra (see Sec- 
tion 0. However, before presenting our results regarding CP Tuc 
modelling, some general results of CYCLOPS should be discussed. 

The process of X-ray emission in polars is isotropic, and the 
emitting region is optically thin; therefore, unlike the optical emis- 
sion, the emitted flux does not change based on the viewing an- 
gle. The orbital variations of the X-ray spectra are explained by the 
temporary obscuration of the post-shock region. This obscuration 
might be due to the accretion column or the white dwarf. In the 
first case, the region that intercepts the sight line is cooler than the 
post-shock region; therefore, it acts as an absorber. This absorp- 
tion is most likely the origin of the broad dip observed in certain 
polars. In the second case, the flux is completely blocked (i.e., the 
flux from the portion of the post-shock region that is located be- 
hind the limb of the white dwarf is not observed). The geometry 
of the mass flux and the orbital inclination completely define both 
cases. Hence, CYCLOPS is a suitable tool for studying these effects 
because it incorporates a 3D representation of the system. 

We chose three sets of parameters to illustrate the two config- 
urations that can produce spectral orbital variations (see TablefTJl. 
To simplify the discussion, we refer to them as absorption, self- 
eclipse I, and self-eclipse II. Figure [T] shows the system in four or- 
bital phases for each model: 0.0, 0.3, 0.5, and 0.8. Figure |2] (left) 
shows the radial variations in density and temperature in the mod- 
els. Figure |2] (right) shows the longitudinal variations in the model 
self-eclipse II. These models are discussed in the following sec- 
tions. 



|http : //heasarc .gsf c ■ nasa.gov/docs/heasarc/caldb/ciocs/memos/cal_ 
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Figure 1. View of the emitting region on the white-dwarf surface for the models absorption (top), self-eclipse I (middle) and self-eclipse II (bottom) - see 
Table[T]- in the orbital phases: 0.0, 0.3, 0.5, and 0.8. Only the walls of the post-shock region (in blue) are represented in this figure. The curved red line is a 
magnetic field line in the accretion column and the green radial line is the magnetic axis. 
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Figure 2. Temperature and density profiles of the models whose parameters are presented in Table \T\ On the left, it is shown the radial profiles used in the 
three models. On the right, the tangential profiles used in the self-eclipse II. 



3.1 Absorption 

Figure |3] (top) shows the phase-resolved spectra for the absorption 
model (see Table [T] and Figure [T] top panel). The spectra tend to 
change along the orbital phase due to the variable absorption that is 
caused by the pre-shock material in the sight line. The absorption 



varies because (i) the quantity of absorbing material varies, and (ii) 
the portion of the obscured emitting region varies. The absorption 
is at its maximum when the region directly points to the observer, 
and most of the emitting region is observed through the pre-shock 
region. This scenario occurs during Phase (see Figures[T] left, and 
Figure[3j magenta short-dashed line). A half- phase later, the region 



Silva et al. 2012 



Absorption 



Table 1. Parameters of the generic models. 



Parameters 



Absorption Self eclipse I Self eclipse II 



i, deg 


33 


^,deg 


18 


^long,deg 


60 


Afl 


0.20 


h,RwD 


0.22 


h 


0.50 


Bpole, MG 


20 


Biat.deg 


74 


Blong, deg 


90 


Tmax,keY 


20 


Nmax, cm-^ (log) 


14.5 


A 


0.26 


B 







0.21 
0. 



58 

68 

60 

0.20 

0.22 

0.50 

20 

45 

360 

20 

14.5 





0.48 



83 

48 

80 

0.10 

0.22 

0.90 

20 

90 

90 

20 

14.5 





0.48 



Af„ 



: mean column density of the pre-shock region (10 cm ). 



points farther from the observer (top panel, Phase 0.5 in FigurelTJ, 
and the pre-shock region is no longer in the observers sight line; 
consequently, the spectrum does not show absorption (green solid 
line in the top panel of Figure |3}. In addition to this geometrical 
effect, the photoabsorption cross-section varies as a function of the 
wavelength. This photoabsorption is larger at lower frequencies, 
which causes a strong change in the continuum emission with re- 
gard to the orbital phase. In particular, the spectra did not show 
difference for energies larger than 3 keV. 
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3.2 Self eclipse 

In the case of self-eclipses, the occulted fraction of the post-shock 
region varied along the orbital cycle. This effect was particularly 
noticeable when the emitting region covered a large area of the WD 
surface. If the temperature changes along the emitting region, the 
partial occultation generates spectral changes. We studied two sit- 
uations that concerned the density and temperature profiles in the 
post-shock region. 

• self-eclipse I - only radial variation; 

• self-eclipse II - only longitudinal variation. 

In the case of self-eclipse I, the WD partially occulted the re- 
gion (see Figure[T] middle panel, Phase 0.5). Consequently, the flux 
decreased for some orbital phases, but significant spectral variation 
was not found (see Figure [3] middle panel). The emission from 
the lower layers (i.e., regions closer to the WD surface) might ex- 
plain the absence of spectral variation, even in the phases when the 
largest occultation of the post-shock region occurred. These regions 
most likely dominated the emission because they are denser that the 
regions closer to the shock front. 

The above result raises the question of whether noticeable 
spectral variations can be obtained when we consider longitudi- 
nal variations in temperature and density. The self-eclipse II model 
(see Table[T]and the lower panel of Figure|3) provides such an ex- 
ample. Using fi — 0.9, the maximum density and temperature of 
the region are located near one of the extremities of the emitting 
region and have values twice as large as the other extremity (Fig- 
ure |2] right). In this model, the temperature and density decrease 
outward in a longitudinal direction, whereas density decreases as 
temperature increases in the radial direction. During Phase 0.5 (see 



Self Eclipse I 




Figure 3. X-ray spectra of the absorption (top), self-eclipse I model (mid- 
dle) and self-eclipse II model (bottom) for 15 equally spaced phases. Some 
phases ai'e overlapped because there is no difference. The spectra ai'e not 
corrected by the instrumental effects discussed in Section [241 



Figure[T] bottom), the hotter and denser region is occulted, whereas 
regions with smaller temperatures and densities remain visible. Fig- 
ure[3](lower panel) shows the spectra along all the cycles. The flux 
levels are different and small changes in the spectral index are also 
observed. Therefore, longitudinal variations in density and temper- 
ature can generate small changes in the shape of spectra. 



3.3 Discussion 

The previous sections described an improvement to CYCLOPS that 
models the X-ray emissions of mCVs. The 3D representation of the 
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Figure 4. (Up) Spectra in arbitrary units for the Absorption model con- 
volved with the ARF and RMF from ASCA satellite. (Bottom) The same for 
the self-eclipse II case. 



accretion column provides a physically consistent description of the 
absorption orbital variation observed in these systems. This new 
representation contrasts with previous models, which relied on ad- 
hoc prescriptions of the orbital energy dependence of absorption. 

Regarding the X-ray orbital variation observed in polars, we 
showed that pre-shock absorption primarily affects the low energy 
region of the spectra, whereas self-eclipse similarly reduces the flux 
of all energies. This difference can be used to distinguish the origin 
of orbital variation in polars. In particular, this absorption might 
better explain the changes in the spectral index of the X-ray spectra 
of polars. 

However, the sensitivity of X-ray instruments might conceal 
the differences described above. To determine whether this hypoth- 
esis is true, Figure |4] shows the generic models convolved with 
ASCA GIS response. This instrument has a low response for en- 
ergies smaller than 1 keV. In this energy range, the differences be- 
tween the absorption and self-eclipse models are larger, but due to 
the instrument response these differences are severely attenuated; 
however, the two models can be distinguished. First, the wave- 
length of the maximum flux in the absorption model moves (in the 
1-2 keV range), whereas this movement is not present in the self- 
eclipse model. Second, the flux does not change for energies larger 
than 3 keV in the absorption model. 



4 APPLICATION: MODELLING CP TUC 

This section describes the simultaneous modelling of the CP Tuc 
optical and X-ray data using CYCLOPS. Section [43] presents the 
previous results concerning CP Tuc that constrain this modelling. 
This section also summarises the previous modelling of this polar. 
Section |42l presents the CP Tuc dataset used in the modelling. The 
application of CYCLOPS to the CP Tuc data is presented in Sec- 
tionl4.3land is discussed in Sectionl4!4l 



4.1 Introduction on CP Tbc 

CP Tu c (AX J23 15- 592) was discovered by the Japanese satellite 
ASCA ( lMisakietaI] |T995). After that report, Thoinas & Reinsc^ 
( Il995h found its optical counterpart and spectroscopically con- 
firmed its classification. A deeper analysis of the optical data re- 
vealed emission lines that contained the two components usually 
found in polars: one narrow line associated with the illuminated 
secondary and one bro ad line originating from th e accretion column 
near the white dwarf (Thomas & Reinschlll99a) . The zero point of 
the spectroscopic ephemeris was defined at the blue-to-red cross- 
ing of the narrow component (i.e., the inferior conjunction of the 
secondary). The broad component presented a maximum blue shift 
near the 0.9 orbital phase, which indicates that this region has the 
smallest angle with regard to the line of sight during this phase. 
iThomas & ReinschI ( Il996p modelled the I b and light curve using 
the cyclotron models of Icharmiugami ( 119921) and found a system 
inclination of 40° and a magnetic field direction in the emitting re- 
gion inclined approximately 30° to the line of sight. They estimated 
the cyclotron spectrum of CP Tuc using the difference between 
the spectra obtained during bright and faint orbital phases. Using 
a temperature of 1 7 keV (as estimated using X-ray data analysis of 
iMisaki et alj ( 119961) - see below), they found that the magnetic field 
was smaller than 17 MG, and the cyclotron emission was optically 
thin in optical wavelengths. 

iMisaki et al.l ( 11996 ) used ASCA data to constmct light curves 
at three X-ray energy intervals. The light curve modulations are 
energy dependent: 87 per cent in 0.7-2.3 keV; 57 per cent in 2.3- 
6 keV; and 14 per cent in 6-10 keV. This type of modulation is 
typical of IPs and usually explained by variable photoabsorption 
along the orbital cycle. The minimum observe d in X-ray co i ncides 
with the inferior conjunction of the secondary. IMisaki et al.l <I996h 
also presented the spectra of CP Tuc combined in two orbital phase 
ranges and fit them using a partial covering model. They suggested 
that the photoabsorption produced by the accretion column modu- 
lates the X-ray light curves. 

iRamsav et alj ltl999) obtained additional X-ray observations, 
optical photometry, and polarimetry for CP Tuc. The absence of 
periodicity other than the orbital period and the presence of white 
light circular polarisation, which reached 10% and was modulated 
by the orbital phase, confiiTned the polar classification; however, 
the energy dependent modulation observed in X-ray light curves 
suggested an IP classification. Significant photometric modulation 
was not found in the B and V bands, but it reached 1.5 (2.0) mag in 
the R (I) band. These findings indicate that th e system has a weak 
magne tic field, as previously suggested by iThomas & ReinschI 
( Il996l) . 

Fixing the fie ld at the magn etic pole at 15 MG and the tem- 
perature at 17 keV. lRamsav et al . ( 1999) modelled the polarimetric 
data and found that the inclination, i, should be larger than 20°. Al- 
though they did not find a unique model that best fit the data, they 
discussed a model with i = 42° and (3 = 50° . The optical light curve 
was modelled considering the self-eclipse of an extended region. 
Moreover, to explain the orbital energy dependence of the X-ray 
light curves, the authors suggested a temperature distribution along 
the longitudinal direction of the emitting region. 

ISeuermann et al.l ( 120071) presented the Zeeman tomography of 
CP Tuc. These observations were conducted in June 2000 when 
the system was in a low state (V=19 m a g). Fixing the parameters 
at those found by iRamsav et al.l (19991), iBeuermann et al.l ( l2007h 
found two equally possible magnetic field configurations: bipolar 
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or multi-polar. In both situations, the magnetic field is low (approx- 
imately 10 MG), and both models reproduce the data equally well. 



4.2 Observational data of CP Tuc 

The optical polarimetric observations and the reduction of CP Tuc 
are presented in Section 14.2.1 1 Our extraction of the spectra using 
the X-ray data obtained bv .Misaki et alj ( ll99q) is described in Sec- 
tion l4.2.2l 



4.2.1 Optical observations 

We obtained the CP Tuc optical data using the Perkin-Elmer 1.6- 
m telescope at the Obser-vatorio do Pico dos Dias (OPD) oper- 
ated by the Laboratorio Nacional de Astroffsica (LNA), using a 
CCD camera modified by the polarimetric module, as described in 
iMagalhaes et alj (199g). The sensor was an EEV front-illuminated 
CCD. Table[2]presents a short description of these observations. 

The redu c tion fo llowed the standard procedures using 
IRAI0 JTodvl Il986l Il993h. The polarisation was calcu- 
lated based on Magalhae s. Benedetti & Roland 11984]) and 
iRodrigues. Cies linski & Steineij d 1998b using the package PC- 
CDPACK (Per eyral 120001) and a set of IRAF routines developed 
by our groupQ The correction to the equatorial reference system 
was performed using standard stars. No instrumental polarisation 
correction was necessary. Our CP Tuc circular polarisation mea- 
surements all had the same sign. However, we were unable to 
calibrate this sign (whether positive or neg ative). We adopted th e 
negative sign based on the measurements of lRamsav et al.l ( 119991) . 

Every linear polarisation measurement has a positive bias: 
The measured polarisation value is greater than the true polarisa- 
tion value (Simmons & Stewart 1985). In particular, measurements 
with P/ap < 1.4 provide only the upper limits of the true polar- 
isation value . The p olarisation of CP Tuc was corrected following 
IVaillancouri ( l2006h . 

The ordinary and extraordinary counts of the polarimetric data 
were summed to obtain light curves. The differential photometry of 
CP Tuc was performed using the star USNO BI.O 0308-0806694 
(R2 = 14.4 mag; I = 13.57 mag) as a reference. The conversion 
between R2 andLandolt's Re for this object indicated a difference 
of 0.03 mag (Kidge? 2003). 

Figure \5\ and |6] present the reduced polarimetric data with 
no binning. From top to bottom, these figures depict magnitude, 
circular polarisation, linear polarisation, and position angle. To 
phase fold our data we used the CP Tuc ephemeris obtained by 
lRamsavetai]h999l) . 

For the CP Tu c modelling (Section [4.3) , we also considered 
iRamsav et al.1 ( Il999l )'s B and V photometry obtained at the SAAO 
I.O-m telescope. 




Figure 5. Optical light curves and polarisation curves of CP Tuc in Re 
band. From top to bottom, magnitude, circular polarisation (V in percent), 
linear polarisation (P in percent), and position angle (6 in degrees). 




Figure 6. Optical light curves and polarisation curves of CP Tuc in Ic band. 
From top to bottom, magnitude, circular polarisation (V in percent), linear 
polarisation (P in percent), and position angle (6 in degrees). 



GIS3 (Table|2l(. CP Tuc was outside the field of view for the Solid- 
iMisakiet al] ( Il996h presented the dataset obtained using the two state Imaging Spectrometer (SIS) instrument. The response matri- 

Gas Imaging Spectrometers (GIS) on board of ASCA: GIS2 and ces were obtained from the ASCA homepage. 

The spectra were extracted at two orbital phase intervals: 0.2- 

0.8, when the system is brighter, and 0.85-0.15, when the system is 
3 IRAF is distributed by National Optical Astronomy Observatories, wliicli fainter. The spectra obtained from GIS2 and GIS3 were combined 

is operated by the Association of Universities for Research in Astronomy, "sing the FTOOL addascaspec and are shown in Figure |7] They 

Inc., under contract with the National Science Foundation do not show significant differences for energies higher than 6 keV, 

* |http: //www .das . inpe .brZ-claudia. rodrigues/polarimetri ghiflaJlieanxffSiJAlblcoflil iron line is seen at approximately 6.4 keV. 
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Table 2. Optical and X-ray data of CP Tuc. 



Date 



Instrament 



Filter T,, 



Duration 



Ref. 



1997 Aug 29 
1997 Aug 30 
1997 Aug 31 
1995 Nov 24/25 
1995 Nov 25 
1995 Nov 02 



Polarimeter + CCD106 + A/4 
Polarimeter + CCD106 + A/4 
Polarimeter + CCD106 + A/2 

Tek8 CCD 

Tek8 CCD 
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This work 
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Figure 7. Phase-resolved spectra of CP Tuc obtained using ASCA data. 
Black: the orbital-phase interval of 0.2-0.8. Red: the orbital-phase interval 
of 0.85-0.15. 



4.2.3 CP Tuc brightness state 

The Re and Ic light curves of lRamsav et alJI 19991 (see their Fig. 1) 
were obtained during the high state and are consistent with our re- 
sults (Figures |6] and |5). The X-ray observat i ons o f Section 14.2.21 
were conducted 23 days before lRamsav et al.l ( ll999l) 's observations 
and were also obtained during the high state. Therefore, we apphed 
our model to high state X-ray and optical data simultaneously. 

We were unable to use the XMM-Newton observations of 
CP Tuc because they were collected during the low state. Accord- 
ing to the AAV SO observation^ the system began in the present 
low state at the end of 2000. In fact, the X-ray observations ob- 
tained at the XMM-Newton obser vatory in 2001 show the system 
in a low state ( Ramsay et alj|2004h . We analysed the data obtained 
from the XMM-Newton observatory in 2004 and 2005, when the 
system was also in a low state. 



4.3 CYCLOPS modelling of CP 'Ric 

As discussed in Section |4~n there are tw o proposed geomet rical 
scenarios for CP Tuc : an absorption mod el (IMisaki et al.l 19961) and 
a self-eclipse model dRamsav et al.l2000l) . The CYCLOPS code han- 
dles both configurations (as discussed in Section[3}, and it is used 
here to unravel the geometrical configuration of the magnetic ac- 
cretion column of CP Tuc. In a preliminary study, we attempted to 
model our polarimetric data and found acceptable solutions using 



|http : //www.aavso .org/l 



any of the CP Tuc proposed geometries dRodrigues et al.l201 ID . We 
called these models Absl and SEl, and their parameters are repro- 
duced in Table[3] 

In the present modelling, we added the light curves in the B 
and V bands from Ramsay et al. (2000) as well as two X-ray spectra 
obtained from Misaki et al.l ([1996) 's data: one in the bright orbital 
phases (0.2-0.8) and the other in the faint orbital phases (0.85- 
0.15). The X-ray reduction is described in Section |4!2l The region 
at approximately 6-7 keV, where a strong iron line is located, was 
not considered in the current fittings. 

We were not interested in fitting the absolute value of flux, so 
we did not use the distance to CP Tuc as an input parameter (eq. 
[8)1. However, we used the same normalisation fcyc across the four 
optical bands (eg. Hit. This method guaranteed that we adjusted the 
total and polarised flux spectral dependence in the optical range. 

The free parameters of CYCLOPS are discussed in Section|2] 
We adopted add equals to 10'^. Due to the low sensitive of ASCA 
GIS in low energy, the spectrum is not sensitive to the B value; 
hence, this parameter was not considered in the CP Tuc fitting. 
We used only one frequency to represent each optical band in the 
search: fi=6.18xl0" Hz, V=5.45xl0^* Hz, 7?c=4.49xl0" Hz 
and /c=3.8x 10^"* Hz. Using many frequencies to represent a band 
does not introduce important differences. In fact, the plotted mod- 
els were calculated considering six frequencies in each band and 
the filter transmission. The X-ray spectra were calculated using 40 
frequencies distributed in the 0.4 to 10 keV interval. 

We determined the best models using two approaches: (i) 
probi ng the optical polarisation models from iRodrigues et al.l 
( I2OII1) , Absl and SEl, and (ii) performing "blind" PIKAIA trials in 
a larger region of the parameter space. The best fit models are Abs2 
and Abs3 (see Figures [8ll9] and [Tol and Table |3}. The Absl, Abs2 
and Abs3 models have absorption geometries, and the SEl model 
presents a self-eclipse of the emitting region. The comparison be- 
tween Figures |4] and [7] suggested that self-eclipse models were not 
appropriate for CP Tuc X-ray data. In fact, the SEl model did not 
reproduce the CP Tuc X-ray spectra (Fig.[TO] dot-dashed line). Nev- 
ertheless, we included the self-eclipse domain in the search to ver- 
ify the robustness of the results of Section [33] But no self-eclipse 
model better than SEl was found. Thus, we discarded self-eclipse 
models for CP Tuc. We henceforth discuss Abs2 and Abs3, which 
have similar x^ values but different geometrical configurations. 

4.4 Modelling discussion 

We begin the discussion of our CP Tuc models by comparing 
certain parameters with previous estimates. Zeeman tomography 
indicated that Bpoie = 19.8 MG (in the bipolar case) and that 
the most frequent magn etic field in the photosphere was 10 MG 
JBeuermann et al-luOOTD . This technique probes the magnetic field 
on the entire WD surface. We found that Bpoie = 6 MG and 7.6 
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Figure 8. Optical light curves of CP Tuc combined in 40 phase bins. From 
top to bottom: B, V, Re and Ic bands. The lines indicate the models Abs2 
(red, dashed l i ne) an d Abs3 (blue, solid hne).The B and V data are from 
lRamsavetai]il999h . 



MG for Abs2 and Abs3, respectively. Table |3] shows the range of 
magnetic fields in the post-shock region (Breg) for each model. 
For instance, Abs3 has a range of 3.8-7.4 MG. Therefore, the mag- 
netic field found using CYCLOPS was smaller that that found by 
ISeuermann et al.l 1 120070 . The CYCLOPS emitting region reached 
a height of 0.14 Rwd, which corresponds to a decrease greater 
than 30% in the magnetic field relative to the value in the pho- 
tosphere. Therefore, the radial extension of the emitting region 
can explain the difference. Zeeman tomography indicated that the 
CP Tu c magnetic field c onfiguration is an offset dipole or a mul- 
tipole JBeuermann et al.ll2007!) . As the CYCLOPS emitting region 
covers a small fraction of the WD surface, our technique is not sen- 
sitive to distinct large-scale magnetic field configurations. 

Misaki et al. ( 1996) estimated a temperature of 17 keV for the 
post-shock region. The maximum temperature of Abs2 and Abs3 
were 60 and 89 keV, respectively. However, the mean temperature 
weig hted by the square d ensity, Tpond, is a more proper comparison 
with lMisaki et al.l y_996)'s single-temperature model. It is 10 keV 
for Abs2 and 15 keV for Abs3. Thus, although we used Tmax to 
define the temperature profile, the hottest regions did not dominate 
the X-ray emission. Given that all the gravitational energy is con- 
verted to X-ray radiation in the post-shock region, we estimated 
the lower limit for the WD mass in CP Tuc using the WD mass- 
radius relationship from iNauenbera (Il972h . We found 0.38 and 
0.47 Mq for Abs2 and Abs3, respectively. The secondary mass was 
0.08 M(T ) from the o rbital period versus the secondary mass rela- 
tion (see lKnigge. Ba raffe & Patterson (2011)). The resulting mass 
ratios were 0.21 and 0. 17 for Abs2 an d Abs3, respectively. 

Optical spectroscopy fo CP Tuc JThomas & Reinschll 19961) is 
one observational constraint that can help to distinguish Abs2 from 
Abs3. The radial velocity of the broad component of the optical 
emission lines displays sinusoidal behaviour with a maximum red- 
shift during Phase 0.9. Hence, the base of the accretion column 
must be observed from the smallest angle at Phase 0.9. Figure [TT] 
shows the system configuration for Abs2 and Abs3 for Phases 0.2, 
0.4, 0.6, 0.8 and 0.9. The curved red line threaded through the emit- 




Figure 9. Polarimetric curves of CP Tuc combined in 40 phase bins for 
Re (top panel) e Ic (bottom panel) bands. The lines indicate the models 
Abs2 (red, dashed line) and Abs3 (blue, solid line). From top to bottom, 
circular polarised flux (Fy), linear polarised flux (Fp) and angle of the 
linear polarisation (6). 



ting region is a magnetic field line, which aids the visualisation of 
the flow direction in the emitting region. The right panels show 
the entire accretion column and illustrate the pre-shock region ge- 
ometry. The system configuration of Abs3 has matched the spec- 
troscopy. We observe the region from the top in Phase 0.9. Thus, 
the post-shock emission has a maximum redshift. The region is 
pointing away from the observer and has the maximum blueshift 
at approximately Phase 0.4. In Abs2, the field line points to the ob- 
server at approximately Phase 0.2 (maximum redshift) and is in the 
opposite direction during Phase 0.8 (maximum blueshift), which 
has not matched the spectroscopy. 

Despite the better agreement between Abs3 and the spec- 
troscopy, we should take this result with caution because it is qual- 
itative. Moreover, the model did not account for certain effects that 
might improve the fit. The Abs2 model presents the best fit of the 
X-ray spectra (see Figure [lO] long-dashed line) and an acceptable 
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Figure 11. Left: View of the emitting region of the models Abs2 (top) and Abs3 (bottom) (Table[T) in five orbital phases: 0.2, 0.4, 0.6 and 0.8 and 0.9. Only 
the laterals (blue) of the post-shock region are represented in this figure. Right: View of the accretion column (orange) for models Abs2 in phase 0.6 (top) and 
Abs3 in phase 0.4 (bottom). In all figures, the curved red line near is a magnetic field line in the accretion column and the green radial line is the magnetic axis. 



Table 3. Parameters of CP Tuc models. 



Parameters 


Abs 1 


SEl 


Abs 2 


Abs 3 


i,deg 


29 


22 


17 


33 


/3,deg 


24 


42 


22 


36 


^long,deg 


18 


13 


16 


32 


Afl 


0.32 


0.52 


0.12 


0.65 


h, RwD 


0.10 


0.22 


0.05 


0.14 


fi 


0.5 


0.5 


0.5 


0.5 


Bp.ie, MG 


6.0 


6.1 


6.0 


7.6 


fi^es.MG 


4.2-5.8 


2.7-5.9 


4.1-5.3 


3.8-7.4 


Biat,deg 


85 


35 


48 


40 


Blong, deg 


41 


336 


63 


337 


imax-> KCV 


66 


102 


60 


89 


T2 ,,keV 

pond 
Nmax,CV[l-^ (log) 


11.0 
15.7 


18.4 
15.0 


10.0 
16.4 


15.0 
13.8 


A^ 


0.26 





0.28 


0.99 


Koi 


3.9 


- 


3.0 


17.0 


Ophase 


-0.002 


-0.095 


0.095 


-0.12 


X^ total 


0.09 


0.22 


0.06 


0.06 


/"P(B),mJy 


12 


11 


12 


12 


/"P(y),mJy 


10 


9 


11 


9 


r^ (Re), m]y 


11 


11 


14 


12 


rnic),miy 


14 


14 


17 


13 



Brcg'- magnetic field range in the post-shock region. 
^ Tpond'- mean temperature weighted using the square density, see Section 
|4.4| for details. 

^ A: the fraction of the maximum possible electron density in the attenua- 
tion region. 
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mean column density of the pre-shock region (10 cm ), 
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0.001 
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Figure 10. X-ray phase resolved spectra in arbitrary units of CP Tuc for 
phase range 0.85-0.15 (top panel) and 0.2-0.8 (lower panel). The lines rep- 
resent the models. 



fit of the optical data (see Figures [S] and |9l long-dashed line). Its 
linear polarisation and cyclotron flux modulation were lower than 
observed. A better agreement between the optical data and Abs2 
might be achieved by including an ellipsoidal variation component 
in the model, because this emission has the proper phasing and 
colour. The elipsoidal variation was not accounted for in the model. 
The Abs3 model (solid lines in Figures [8] l9l and [TOt provides 
the best fit to the optical data, with linear polarisation and flux mod- 
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ulation consistent with the observations. It also presented an accept- 
able fit to the high-energy X-ray spectrum, but it overestimated the 
absorption at low energies. This excess absorption might be due to 
the adopted approximations of the pre-shock material, which was 
cold, homogeneous, and had solar abundance. A warmer absorber 
with smaller abundance might have less photoabsorption. A decre- 
ment in the density of regions farther from the WD might also ac- 
count for the smaller net absorption. 

Importantly, CYCLOPS only includes the bremsstrahlung con- 
tinuum emission in the X-ray range, which is the dominant process. 
However, other processes might als o contribute includ ing Compton 
scattering in the post-shock region (ISuleimanov et a l. 2008), heat- 
ing of the WD surface, and emission lines. 

Taking into account all of the above considerations, small dif- 
ferences between the observed data and the model were expected. 
On the other hand, the quantity and quality of the dataset might 
be not enough to justify the inclusion of these "second-order" pro- 
cesses in the models. However, we argue that these data are enough 
to restrict the correct viewing aspect of this system, which has 
phase-resolved absorption in X-rays. 



5 CONCLUSIONS 

We extended the CYCLOPS code to study the X-ray and optical con- 
tinua emission of polars. We included bremsstrahlung X-ray emis- 
sion, the dominant radiative process in the post-shock region, and 
photo-absoiption by the pre-shock accretion column. Cyclotron 
emission and self-eclipse by the white dwarf were already avail- 
able in the previous version of the code. The new version allowed 
us to simultaneously model X-ray and optical data. 

The calculated spectra accounted for the temperature and den- 
sity structure in the post-shock region in both radial and tangential 
directions. Although the X-ray emission of polars is usually op- 
tically thin, our radiative transfer solution adequately treated op- 
tically thick and thin regions. The magnetic field lines in the 3D 
space defined the accretion column geometry. This approach re- 
sulted in a consistent representation of the geometrical effects of 
the white dwarf occultation and the pre-shock region absorption. 
Moreover, the code provided phase-resolved X-ray spectra. 

We showed that phase-resolved X-ray spectra have distinct 
signatures in cases of self-eclipse and absorption. Absorption 
causes a clear variation of the spectral index along the orbital cy- 
cle, whereas these differences are barely distinguished in the self- 
eclipse case, even when we consider the inhomogeneous distribu- 
tion of the temperature and density in the emitting region. 

We chose CP Tuc to explore the differences between self- 
eclipse and absorption in the X-ray spectra of polars because it 
is a polar with two diverging explanations for its X-ray data, and 
each one is consistent with one of the above scenarios. We mod- 
elled CP Tuc using the X-ray and optical data from the literature as 
well as new optical Re and /c polarimetry and photometry data. 

Our CP Tuc models indicated a single accretion region ob- 
ject in which cyclotron beaming generated the optical modulation 
and photo -absorption generate d the X-ray modulation, in agree- 
ment with lMisaki et alj JI996l) . We could not create a model that 
described the X-ray spectra using the self-eclipse scenario. The 
magnetic field intensity of our best fit tings was consistent wit h pre- 
vious Zeeman tomography estimates dBeuermann et aLluOOTT) . The 
temperature stmct ure matched the sin gle-temperature model of the 
X-ray continuum JMisaki et alj|l996h . One of the best fittings had 
the correct phasing to account for the radial velocity curve of the 



broad component of the optical emission lines JThomas & ReinschI 
[1993). 

Future improvements to the code are a more physically consis- 
tent description of the pre-shock and post-shock density and tem- 
perature distributions and inclusion of other emitting process. With 
regard to CP Tuc, a more detailed model requires more quality data 
in high-energy range. 
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